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IV 


ABSTRACT 


We  develop  several  models  for  the  spread  of  two  opposing  ideologies  in  a  elosed 
population  based  on  epidemie  models.  Based  on  different  interaction  rules  for  the 
mechanism  of  how  the  ideologies  are  spread,  we  study  the  results  of  deterministic  and 
stochastic  models.  The  goal  of  our  work  is  to  provide  a  tractable  analytical  framework  for 
each  hypothesized  set  of  social  interactions,  and  to  analyze  the  effect  of  different  initial 
conditions  on  the  proportion  of  the  population  affiliated  with  each  ideology  after  a  large 
time  interval.  The  models  developed  herein  are  designed  to  give  decision  makers  an 
insight  into  a  complex  process. 
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EXECUTIVE  SUMMARY 


One  of  the  objectives  listed  in  the  2003  National  Strategy  for  Combating 
Terrorism  is  to  win  the  “War  of  Ideas”.  i  Our  work  seeks  to  place  an  analytic  framework 
around  this  war.  Our  goal  is  to  create  a  methodology  for  considering  alternatives  and 
some  concrete  metrics  with  which  to  compare  courses  of  action.  Our  fundamental 
assumption  is  that  one-to-one  (i.e.  interpersonal)  communication  is  the  most  important  in 
spreading  ideas. 

Our  tools  are  deterministic  and  stochastic  models  originally  developed  for 
infectious  diseases  and  rumor  propagation.  The  fundamental  idea  is  that  when  two 
people  in  the  population  connect  (either  through  direct  contact  or  phone,  email,  etc),  the 
ideology  may  be  spread.  These  models  are  similar  to  traditional  epidemic  models.  Many 
extensions  to  the  idea  of  the  spread  of  ideology  as  disease  are  possible  and  some  are 
explored  in  our  work. 

We  extend  previous  work  by  placing  ideology  in  a  greater  social  context.  We 
introduce  two,  diametrically  opposed  ideas  and  model  their  flow.  We  refer  to  the 
proponents  of  these  ideas  as  the  supporters  and  contrarians.  We  consider  the  case  where 
both  the  supporters  and  contrarians  openly  vie  for  a  greater  share  of  support  from  the 
public.  This  case  is  similar  to  a  political  campaign  in  the  United  States  without  the 
influence  of  media.  We  also  consider  the  case  where  the  supporters  are  able  to  openly 
propagate  their  message,  but  the  contrarians  only  interact  when  supporters  try  to  convert 
them.  We  believe  this  is  the  case  where  there  is  a  small  but  dedicated  extremist  sub¬ 
population.  We  find  that  under  the  model  assumptions  a  relatively  small  number  of 
contrarians  are  required  to  overcome  a  large  increase  in  the  supporters 


1  National  Strategy  for  Combating  Terrorism,  23. 
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I.  INTRODUCTION 


From  the  beginning,  the  War  on  Terror  has  been  both  a  battle  of  arms  and 
a  battle  of  ideas  -  a  fight  against  the  terrorists  and  against  their  (ideology). 

In  the  short  run,  the  fight  involves  using  military  foree  and  other 
instruments  of  national  power. . .  In  the  long  run,  winning  the  war  on  terror 
means  winning  the  battle  of  ideas,  for  it  is  ideas  that  can  turn  the 
disenchanted  into  terrorists. 

—  National  Security  Strategy  of  the  United  States,  March  2006 

It  is  ideological  belief,  reinforced  by  propaganda  operations,  that 
convinces  recruits  and  supporters  that  their  actions  are  morally  justifiable. 

All  the  instruments  of  national  power  play  a  role  in  undermining 
ideological  support...  Targeting  ideology  includes  amplifying  the  voices 
of  those  who  promote  alternative  ideas 

—  National  Military  Strategic  Plan  for  the  War  on  Terrorism, 
February  2006 

One  of  the  objectives  listed  in  the  National  Strategy  for  Combating  Terrorism 
(February  2003)  is  to  “Win  the  War  of  Ideas”.  Specifically: 

Together  with  the  international  community,  we  will  wage  a  war  of  ideas  to 
make  clear  that  all  acts  of  terrorism  are  illegitimate,  to  ensure  that  the 
conditions  and  ideologies  that  promote  terrorism  do  not  find  fertile 
ground...  and  to  kindle  the  hopes  and  aspirations  of  freedom  of  those  in 
societies  ruled  by  the  sponsors  of  global  terrorism. 

Clearly,  the  United  States,  at  the  highest  level,  values  winning  the  War  of  Ideas, 
and  for  good  reason.  Winning  the  war  of  ideas  is  possibly  the  single  most  important 
action  the  US  can  take  to  further  the  National  Interest.  Winning  the  War  of  Ideas  would 
have  the  immediate  benefit  of  significantly  reducing  the  requirement  for  military  force. 
Winning  the  War  of  Ideas  will  cost  far  less  than  imposing  our  will  by  force. 

To  our  knowledge,  little  work  has  been  done  by  Operations  Research  practitioners 
on  this  important  problem.  At  first  glance,  the  statement  of  this  problem  lends  itself  more 
readily  to  the  social  scientists.  Operations  Research  practitioners  do  in  fact  have 
significant  tools  to  contribute.  The  nature  of  the  current  war  is  new  and  different  than 
past  conflicts  and  requires  all  of  us  to  step  outside  our  traditional  boundaries. 
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Our  examination  of  the  war  of  ideas  takes  plaee  at  the  miero  level,  foeusing 
exclusively  on  one-to-one  communications.  These  may  be  instantiated  either  in  person  or 
via  a  communication  channel  such  as  email  or  telephone.  We  do  not  consider  the  effect 
of  mass  media.  Apart  from  mass  media,  we  believe  that  the  most  important 
communication  mode  in  swaying  a  person’s  ideology  is  interpersonal  communications. 
We  use  a  simple  model  to  describe  a  person’s  ideology;  they  can  either  be  neutral,  for,  or 
against.  The  object  of  their  ideology  is  immaterial  for  our  purposes.  We  view  the  spread 
of  ideology  to  be  an  epidemic  with  basically  two  mutually  exclusive  populations: 
supporters  and  contrarians.  We  extend  classic  models  of  epidemiology  and  rumor 
propagation  to  suit  this  new  purpose. 

Our  goals  in  this  thesis  are  threefold.  First,  we  seek  to  provide  an  analytic, 
descriptive  framework  for  decision  making  in  the  War  of  Ideas.  To  clarify,  we  will 
highlight  what  our  model  is  not:  It  is  neither  a  proscriptive  set  of  equations  to  be  used  in 
a  vacuum,  nor  a  predictive  model  that  projects  population  behavior.  Rather,  it  is  a 
descriptive  tool  to  help  understand  cause-and-effect  relations  in  the  context  of  ‘war  of 
ideas’.  Secondly,  we  provide  coarse-detail,  deterministic  and  stochastic  models  for  the 
propagation  of  ideology,  which  may  be  used  as  a  base  case  for  validating  complex 
simulations.  The  US  Army  is  currently  in  the  process  of  examining  proposals  for 
simulation  models  to  study  the  War  of  Ideas2.  Finally,  our  model  is  mathematically 
interesting  because  it  introduces  competition  in  the  context  of  rumor  propagation.  In 
addition  to  being  interesting  from  a  purely  mathematical  point  of  view,  the  models 
developed  herein  represent  a  step  forward  in  the  study  of  social  interactions  as  stochastic 
processes. 

Chapter  II  provides  a  short  review  of  current  models  for  the  spread  of 
rumor/ideology.  Chapter  III  develops  our  models,  called  NSRL  and  NSRCL  and  explores 
their  behavior  deterministically.  Chapter  IV  explores  the  NSRL  model  probabilistically 
using  both  a  Markov  chain  as  well  as  a  diffusion  approximation.  Chapter  V  presents  our 
conclusion  and  recommends  future  research. 


2  U.S.  Army  Small  Business  Innovative  Research  initiative  A06-094 
http://www.acq.osd.mil/osbp/sbir/solicitations/sbir062/army062.htm 
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II.  LITERATURE  REVIEW 


The  mathematical  theory  of  epidemics  appears  at  present  to  be  developing 
fairly  rapidly.  The  theory  is  only  likely  to  have  valuable  applications  in  so 
far  as  it  is  developed  in  the  context  of  a  proper  understanding  of  the 
epidemiological  realities.  Many  of  the  models...  are  of  necessity 
oversimplified:  nevertheless,  useful  results  are  readily  available. 

-  Norman  T.  J.  Bailey,  1975. 

A,  RUMORS  AS  A  SUBSET  OF  INFECTIOUS  DISEASE  MODELS 

The  analytic  models  of  rumor  (used  interchangeably  with  information) 
propagation  are  closely  linked  to  the  analytic  models  of  infectious  diseases.  Indeed,  one 
of  the  earliest  works  in  this  field  was  motivated  by  the  desire  to  control  ‘psychological 
infections’  and  was  developed  in  the  1920’s  (Daley  and  Gani,  1999).  The  similarity 
between  rumor  propagation  models  and  models  for  infectious  disease  is  that  they  both 
describe  the  behaviors  of  actors;  we  use  the  term  ‘actor(s)’  in  our  development  to  mean 
the  atomic  level  of  social  interaction.  Actors  may  be  thought  of  as  individuals,  but  may 
also  be  thought  of  as  family  units  or  even  entire  villages.  This  definition  of  an  actor  is 
encouraged  both  in  Maki  and  Thompson  (1973)  as  well  as  by  Wasserman  and  Faust 
(1994).  The  development  of  models  for  infectious  disease  focuses  on  these  actors 
interacting  with  each  other  in  accordance  with  some  well-defined  (deterministic  or 
stochastic)  set  of  rules.  In  this  context  it  is  useful  to  think  of  spreading  a  rumor  as 
‘infecting  with  knowledge’.  We  review  two  models  for  the  spread  of  rumors:  The 
Maki-Thompson  [MT]  (Maki  and  Thompson  (1973))  and  Daley-Kendall  [DK]  (Daley 
and  Kendall  (1964))  models. 


B,  MODELS  FOR  INFECTIOUS  DISEASE:  THE  BASIC  FRAMEWORK 

Since  our  development  depends  heavily  on  the  existing  mechanics  of  infectious 
disease  models,  it  is  useful  to  describe  how  infectious  disease  models  work.  Particularly 
useful  subsets  of  infectious  disease  models  are  those  which  partition  the  population  into 
three  classes.  The  following  table  explains  the  so-called  “S-I-R”  (after  Waltman,  1970) 
models.  A  rigorous  development  of  the  basic  mathematical  models  for  infectious 
diseases  is  contained  in  Bailey  (1975). 
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Population 

Definition 

S 

The  Susceptible  population.  These  are  individuals  who  are  not  yet 

infected,  but  may  be  at  some  future  time.  In  rumor  models  they  are 

called  the  ignorant  population. 

I 

The  Infective  population.  These  individuals  are  the  active  carriers  and 

spreaders  of  disease.  In  rumor  models  they  are  called  spreaders. 

R 

The  Removed  population.  These  individuals  have  been  infected  and 

have  since  been  removed  (either  by  death  or  immunity).  In  rumor 

models  they  are  called  stiflers. 

Table  1:  Definitions  of  States 


A  eharacteristic  of  S-I-R  models  is  a  finite  population  of  size  N  +  1,  representing 
N  initial  suseeptibles  and  one  initial  carrier.  Furthermore,  this  population  exists  within  a 
closed  system.  Because  the  future  state  of  the  system  is  dependant  only  upon  the  current 
state,  this  model  exhibits  the  Markov  property. 

This  type  of  model  belongs  to  the  class  of  strictly  evolutionary  processes  because 
of  the  following  property:  once  a  state,  as  defined  by  the  vector  ,  is  visited 

it  may  never  be  visited  again.  In  both  the  [DK]  and  [MT]  models,  which  are  discussed  in 
the  next  section,  the  members  of  the  population  who  make  the  transition  from  S  to  I 
eventually  reach  the  absorbing  state  of  R.  For  both  deterministic  and  stochastic  models, 
at  the  end  of  the  rumor  spread,  the  population  is  partitioned  into  ignorants  (S)  and  stiflers 
(R),  with  no  spreaders  (I)  left.  In  the  stochastic  version  of  this  model,  an  embedded 
Martingale  is  present  as  noted  by  Daley  and  Gani,  (1999).  We  use  the  embedded 
Martingale  in  Chapter  IV  to  study  the  asymptotic  distribution  of  the  process  without 
explicitly  specifying  the  associated  Markov  chain. 

Belen  and  Pearce  (2004)  note  that  formally,  the  models  described  by  the  two 
models  ([DK]  and  [MT])  are  the  same. 
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C.  RUMOR  PROPAGATION 
1.  The  Basic  Model 

The  basic  mechanics  of  rumor  propagation  are  described  in  (Daley  and  Kendall, 
1965)  and  are  as  follows:  The  population  is  partitioned  into  three  mutually  exclusive  and 
collectively  exhaustive  classes  as  in  Table  1.  The  classes  within  the  population  have  the 
following  interpretation: 

S:  Susceptibles  (interchangeably,  the  ignorants).  Those  who  do  not  know  the 
information  but  may  in  the  future. 

I:  The  active  spreaders.  Those  who  know  the  information  and  are  actively 
disseminating  it  to  those  individuals  they  meet 

R:  The  stiflers.  The  individuals  who  know  the  information  and  are  no  longer 
spreading  it. 

The  same  definitions  are  used  in  the  [MT]  model.  A  stochastic  version  of  this 
same  model  is  provided  as  Appendix  B.  The  properties  of  the  [DK]  model  have  been 
explored  by  several  authors,  including  Watson  (1987),  Pittel  (1990),  and  Hayes  (2005). 


The  following  table  describes  the  transitions  between  actor  A  and  actor  B: 


ACTOR “A” 

S 

I 

R 

S 

No  Change 

B->I 

No  Change 

ACTOR “B” 

A:  no  change 

I 

A->I 

A->R 

B->R 

B->R 

R 

No  Change 

A->R 

No  Change 

Table  2:  Daley-Kendall  Transition  rules  (after  Hayes,  2005) 


The  [MT]  model  is  very  similar,  except  that  when  A  and  B  are  both  members  of 
class  I,  only  one  becomes  a  stifler  as  a  result  of  the  interaction. 

Equations  (1)  present  the  deterministic  Daley-Kendall  model  (Daley  and  Gani, 

1999): 
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S'=-SI 

r=si-i^-iR 

R'  =  I^+IR 


(1) 


The  numerical  solution  of  equations  (1  is  presented  as  Figure  (1)  with  a 
population  size  of  1000,  and  1  initial  spreader. 
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Figure  1.  Daley-Kendall  Model 

2,  Stifling  Mechanisms 

Common  experience  tells  us  that  rumors  do  not  generally  spread  indefinitely.  To 
force  the  rumor  to  end  at  some  time,  I  must  be  a  transient  state.  The  mechanism  by 
which  actors  leave  the  class  I  is  called  stifling. 

a.  Models  Which  Lack  Stifling  Mechanisms 

The  simplest  models  lack  stifling.  An  example  of  this  type  of  model  is 
attributed  to  Bailey  (1957)  and  is  described  in  chapter  9.2  of  Bartholomew  (1967).  This 
model  can  be  derived  directly  from  [DK]  by  neglecting  the  class  R.  In  a  deterministic 
model  which  lacks  stifling,  the  number  of  spreaders  increases  as  a  logistic  ^-curve  until 
only  one  ignorant  remains\  At  this  stage  the  rumor  is  considered  to  have  ‘fully 
propagated’.  In  the  stochastic  model,  the  rumor  is  fully  propagated  when  no  ignorants 
remain. 

3  In  the  deterministic  model,  the  number  of  ignorants  asymptotically  reaches  zero.  Therefore,  a 
limiting  argument  on  one  ignorant  remaining  must  be  used. 
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In  the  no-stifling  formulation,  the  population  consists  only  of  S  and  / 
types,  and  because  we  assume  a  known  fixed  population  N,  this  model  collapses  to  a  one 
dimensional  case  because  S{t)  =  N  -  I(t).  This  finite  population  property  is  frequently 
exploited  in  the  subsequent  Chapters  of  this  Thesis  to  simplify  calculations.  Equations 
(2)  and  (3)  contain  the  equation  for  the  deterministic  simple  epidemic  in  continuous  time 
(Daley  and  Gani,  1999)  and  its  solution,  respectively.  Figure  (2)  presents  a  graph  of 
equation  (3) 


dl 

dt 


=  pSI=/3I{N-I) 


(2) 


^(0 


IqN 


(3) 


time 


Figure  2.  Deterministic  non-stifling  sample  path 

b.  Daley-Kendall  Stifling 

Everyday  experience  tells  us  that  not  all  news  and  rumors  will  propagate 
forever.  A  stifling  mechanism  is  necessary  if  the  rumor  is  to  spread  for  a  finite  period  of 
time  (Daley  and  Gani,  (1999)). 

In  Daley  and  Kendall’s  model,  stifling  occurs  when  a  spreader  meets 
another  spreader,  or  when  a  spreader  meets  a  stifier.  In  the  instance  of  the  [DK]  model. 
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when  two  spreaders  meet,  both  beeome  stiflers.  In  the  instanee  of  the  [MT]  model,  both 
only  one  spreader  beeomes  a  stifler. 


D.  THEORETICAL  RESULTS 

An  interesting  result  was  first  demonstrated  in  the  deterministic  version  of  the 
[DK]  model,  later  extended  to  the  stochastic  version  and  finally  (Belen  and  Pearce  2004) 
generalized  to  the  [MT]  model:  Given  one  initial  spreader,  the  proportion  of  ignorants  at 
the  end  of  the  rumor  propagation  converges  to  approximately  20%.  This  result  is 
published  in  Daley  and  Kendall  (1965),  and  is  confirmed  in  Pittel  (1990)  as  well  as  Hayes 
(2005)  and  is  mentioned  by  Zanette  (2005).  A  proof  is  contained  in  the  aforementioned 
papers,  as  well  as  in  Daley  and  Gani  (1999).  Additionally,  as  shown  in  Daley  and  Gani 
(1999),  for  ‘small’  populations,  there  is  a  probability  of  very  limited  spread  due  to  the 
spreaders  encountering  other  spreaders  immediately  at  the  start  of  rumor  propagation.  As 
the  population  tends  towards  infinity,  the  probability  of  this  small  spread  phenomenon 
tends  towards  zero  (Gani  and  Daley,  1999). 

The  small  spread  phenomenon  is  of  great  interest  when  studying  the  stochastic 
behavior  of  these  processes.  We  will  discuss  this  issue  extensively  in  Chapter  IV. 

E.  POPULATION  DYNAMICS 

The  epidemic-style  models  of  rumor  propagation  assumes  a  homogeneous  mixing 
of  the  population  and  results  in  diffusion  according  to  the  Law  of  Mass  Action  (Daley 
and  Gani,  1999).  Other  population  mixing  dynamics  are  possible  and  have  been  explored 
by  several  authors,  c.  f.  Kress  (2005).  In  the  absence  of  any  information  about  the 
interaction  between  the  actors,  homogeneous  mixing  is  an  appropriate  assumption.  With 
limited  information  about  the  sociodynamics  of  the  population,  mechanisms  like  those 
proposed  by  Kress  (2005)  are  appropriate. 

If  more  detailed  information  about  the  underlying  social  organization  exists,  then 
a  different  approach  is  warranted.  Zanette  (2005)  explores  the  propagation  of  rumors  on 
both  ‘small  world  networks’  (deterministic)  as  well  as  ‘dynamic  small  worlds’ 
(stochastic).  Zanette’s  approach  builds  both  on  the  [DK]  model. 
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F.  DETERMINISTIC  VS.  STOCHASTIC  MODELS 

The  model  of  rumor  propagation  may  be  explored  by  either  deterministie  models 
or  stoehastie  models.  The  approaeh  taken  depends  on  the  level  of  analysis  the  researeher 
seeks  to  eonduet  and  also  depends  to  some  extent  on  the  eomplexity  of  the  problem.  A 
deterministie  modeling  approaeh  may  be  taken  either  in  diserete  or  eontinuous  time.  The 
diserete  time  model  is  implemented  as  a  set  of  differenee  equations,  deseribing  system 
behavior  at  time  0,  At ,  2At ,. . . .  The  limit  of  the  diserete  time  model  as  At  approaehes 
zero  is  a  eontinuous  time  model.  The  eontinuous  time  model  is  implemented  as  a  set  of 
differential  equations.  With  either  the  diserete  or  eontinuous  time  deterministie  models, 
there  will  be  only  one  sample  path  for  eaeh  set  of  initial  eonditions  and  parameter  values. 
Given  the  state  veetor  of  the  number  of  individuals  in  the  population  in  eaeh  state  (S,  I, 
and  R)  at  any  time,  the  future  state  veetor  of  the  number  of  individuals  in  eaeh  state  is 
eompletely  speeified. 

Stoehastie  models  take  a  different  approaeh.  The  stoehastie  proeess  in  diserete 
time,  modeled  by  a  diserete -time  Markov  ehain  foeuses  on  the  probabilities  of  transition 
between  states.  In  sharp  eontrast  to  the  deterministie  model,  knowing  the  state  veetor  of 
the  number  of  individuals  in  eaeh  state  (S,  I,  R)  at  any  time  does  not  (unless  the  state 
deseribed  is  absorbing)  speeify  a  single  future  state  of  the  system.  Many  sample  paths 
are  possible.  The  stoehastie  approaeh  will  yield  eonditional  probabilities  of  future  states 
given  the  eurrent  state.  As  we  are  only  studying  the  asymptotie  behavior  of  the  proeess  as 
^  T  00 ,  we  do  not  have  to  make  any  speeifie  statement  about  the  distribution  of  transition 
times. 

Although  generally  more  diffieult  to  implement  and  analyze,  stoehastie  models 
provide  a  riehness  that  deterministie  models  laek.  Stoehastie  models  give  us  not  only  the 
expeeted  outeome  of  the  rumor  proeess,  but  the  distribution  of  outeomes  as  well. 
Deterministie  models  have  the  advantage  of  being  easy  to  formulate  and  implement. 
They  are  also  the  only  ehoiee  in  situations  where  the  state  spaee  that  arises  in  the 
stoehastie  formulation  is  simply  intraetable.  The  following  quote  from  Bailey  (1975) 
shapes,  in  our  opinion,  the  trade-offs  between  purely  deterministie  and  stoehastie  models; 
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The  mathematician’s  first  choice  was  always  to  go  for  stochastic 
modeling,  since  it  was  more  interesting;  and  in  any  case  there  were  plenty 
of  examples,  such  as  the  analysis  of  family  data,  for  which  deterministic 
models  would  be  virtually  meaningless.  [However]  ...  a  somewhat  more 
balanced  judgment  is  possible.  When  numbers  of  susceptibles  and 
infectives  are  both  large  and  mixing  is  reasonably  homogenous,  a 
deterministic  treatment  is  likely  to  be  fairly  satisfactory  as  a  first 
approximation,  though  even  here  it  may  be  difficult  to  account  for  large- 
scale  phenomena. .  A 


4  Bailey,  .  The  Mathematical  Theory  of  Infectious  Diseases  ,  25. Bailey,  Norman  T.J. 
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III.  DETERMINISTIC  MODELS 


A.  MOTIVATION,  PRELIMINARIES  AND  NOTATION 

Deterministic  models  in  the  form  of  systems  of  differential  equations  are  of  great 
importance  in  the  study  of  evolutionary  processes  such  as  ideological  propagation  and 
epidemics.  Several  texts,  notably  Anderson  and  May  (1992),  are  devoted  solely  to  their 
development  and  implementation.  A  deterministic  model  is  one  where,  at  each 
infinitesimal  time  step,  the  underlying  (probabilistic)  process  is  assumed  to  take  on  its 
expected  value,  resulting  in  a  single  path  through  the  state  space. 

This  class  of  models  is  attractive  for  several  reasons.  First,  they  are  easy  to 
implement.  If  the  given  system  does  not  yield  a  closed-form  solution,  we  may  still 
investigate  its  asymptotic  (long  range)  behavior  or  solve  the  system  numerically  using 
one  of  many  available  methods.  In  general,  numerical  solution  is  very  fast. 

Because  deterministic  models  do  not  respect  the  integer  nature  of  populations  as 
stochastic  models  do,  their  qualitative  behavior  is  invariant  to  the  population  size.  This 
feature  makes  them  unsatisfactory  for  small  populations  where  random  fluctuations  may 
be  significant  and  therefore  stochastic  models  should  be  used.  However,  in  Chapter  IV 
we  reference  a  theoretical  result  (Barbour,  1972),  which  shows  that  under  certain 
assumptions  for  large  populations,  the  stochastic  mean  converges  to  the  deterministic 
solution. 

The  models  we  develop  in  the  following  section  will  be  used  in  this  Chapter  and 
Chapter  IV.  In  these  models  we  identify  several  population  classes  that  interact.  For  any 
population  class  A  we  denote  the  initial  population  size  as  Xo,  and  the  long  term  limiting 
size  (if  such  limit  exists)  asX^ .  We  use  X=  X{t)  to  denote  also  the  size  of  the  population 
at  time  t.  In  general  we  consider  four  population  classes: 

•  N:  Neutrals  who  hold  no  firm  ideology  or  opinion  and  who  may  be 
susceptible  to  ideological  influence; 

•  S:  Active  Supporters  who  actively  spread  pro-government  ideas; 
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•  R:  Latent  Supporters  who  support  the  government  but  do  not  aetively  spread 
their  ideas; 

•  C;  Active  Contrarians  who  aetively  spread  anti-government  ideas; 

•  L:  Latent  Contrarians  who  oppose  the  government  but  do  not  aetively  spread 
their  ideas; 

•  P=A^+5'+i?++C+Z;  Total  Population. 

In  Seetion  B  we  develop  two  models  ealled  NSRL  and  NSRCL.  In  Seetion  C  we 
examine  some  speeial  eases  of  these  models  whieh  possess  elosed-form  solutions.  In 
Seetion  D  we  present  numerieal  results. 

B,  THE  NSRL  AND  NSRCL  MODELS 

We  develop  two  models  whieh  eapture  the  dynamies  of  the  spread  of  ideas. 
Although  the  two  models  represent  different  situations,  they  share  a  similar  strueture. 
Both  models  are  posed  here  in  a  deterministic  form  for  closed  populations;  that  is, 
individuals  do  not  enter  or  leave  the  population.  This  assumption  is  realistic  in  remote 
isolated  communities.  A  generalized  case,  allowing  for  migration  in  and  out  of  the 
population  is  a  natural  and  feasible  extension,  but  it  is  one  we  do  not  pursue  here. 

1.  The  AAPL  Model 

The  NSRL  model  considers  the  following  situation:  There  is  a  message,  which  is 
overtly  being  passed  by  the  active  supporters  (A).  There  are  a  certain  number  of  latent 
contrarians  (L)  in  the  population  who  are  antagonistic  to  the  supporters,  but  who  do  not 
actively  spread  their  views.  These  contrarians  are  passive  until  they  meet  an  active 
supporter,  in  which  case  they  will  attempt  to  convert  him  to  their  side.  The  implication  is 
that  individuals  with  a  counter-attitude  do  not  want  to  advertise,  but  when  an  active 
supporter  attempts  to  recruit  a  contrarian,  the  counter-message  by  the  contrarian 
overwhelms  (with  some  probability)  and  converts  the  supporter  to  the  contrarian  way  of 
thinking.  The  NSRL  model  is  asymmetrical  in  the  sense  that  the  supporters  and 
contrarians  sides  behave  differently  (e.g.,  there  are  no  overt  contrarians)  as  opposed  to  a 
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symmetric  model,  in  which  both  sides  have  overt  members.  We  have  the  situation  of  an 
insurgency  or  popular  front  in  mind  in  the  development  of  the  NSRL  model. 

Initially  the  majority  of  the  population  is  in  class  N  (neutrals)  with  some 
proportion  belonging  to  class  L  and  some  proportion  belonging  to  class  S. 

We  assume  that  homogeneous  mixing  takes  place  in  the  population,  meaning  that 
the  probability  of  interaction  between  two  members  is  dependant  on  the  proportions  and 
attributes  of  the  populations  only,  and  not  other  factors  such  as  time  of  day,  location,  etc. 
An  extension  involving  time  variant  mixing  parameters  is  feasible,  but  not  pursued  here. 

An  interaction  occurs  when  one  member  (of  any  class)  meets  another  member 
(also  of  any  class),  and  the  result  of  this  meeting  is  that  at  least  one  of  the  two  members 
may  change  its  class.  Table  3  presents  the  seven  possible  interactions  in  the  NSRL  model. 
The  bold  letters  in  the  table  are  the  row  and  column  indices  (classes).  The  entry 
corresponding  to  the  row/column  pair  describes  the  possible  class  change  that  may  occur. 
For  example,  ‘N->S’  means  that  a  member  of  the  N  population  class  may  leave  its  class 
and  join  the  S  class.  We  will  use  this  notation  for  the  remainder  of  the  chapter. 


N 

5 

R 

L 

N 

— 

N->S 

— 

— 

S 

N->S 

S->R 

S->R 

S->R 

S->L 

R 

— 

S->R 

— 

— 

L 

— 

S->L 

— 

— 

Table  3;  NSRL  Interactions 
Figure  3  is  a  pictorial  representation  of  Table  3. 
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Figures.  NSRL  Interactions 
It  is  evident  that  class  S  interacts  with  every  other  class  (including  itself).  All 
other  classes  only  interact  with  class  S.  Table  3  leads  directly  to  the  following  system  of 
Ordinary  Differential  Equations  (ODEs): 


—=-f3,NS 

dt 

— = J3,NS  -Jd^S^-  J3,SR  -  /3,SL 
dt 

dt 


with  the  initial  conditions: 
S,=r,N 

Ro-0 

L,=r,N 
fs  +/y<l 


(4) 


where  the  are  parameters  describing  the  intensity  of  each  class’  interaction  with  the 
other  classes,  Nq,Sq,Rq,Lq  are  the  initial  sizes  of  each  class,  and  r^,r^  are  the  proportion 
of  the  population  that  is  initially  in  classes  S  and  L,  respectively. 
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Note  that  it  is  possible  to  fix  the  parameter  P  and  appropriately  adjust  the  |3^  s  to 

achieve  the  same  result.  We  will  employ  this  technique  and  use  the  ‘normalized’  system 
with  P=\  in  the  analysis  section. 

We  allow  for  the  possibility  that  members  of  different  populations  may  act  with 
different  intensities.  For  example,  the  active  supporters  {S)  may  aggressively  seek  out 
Neutrals  (TV),  but  avoid  Retired  (if)  supporters.  Note  that  if  =  0  ,  the  model  posed  here 

collapses  to  that  proposed  by  Daley  and  Kendall  (1964).  If  =  0 ,  then  the  initial  state  is 
an  absorbing  state  and  no  change  takes  place. 

2  The  Model 

This  model  has  a  similar  development  to  the  NSRC,  with  one  important 
difference:  it  is  symmetrical.  We  add  the  class  ‘C’  of  active  contrarians  and  now  the  two 
sides  have  symmetric  rules. 

The  situation  described  by  this  model  differs  from  the  asymmetrical  NSRL  model 
in  that  both  messages  are  spread  openly.  This  situation  may  be  one  of  political  change, 
such  as  the  formation  of  a  government  or  an  election.  The  active  supporters  and  active 
contrarians  may  be  thought  of  as  belonging  to  different  political  parties,  voting  for 
different  candidates,  etc.  The  interactions  in  this  model  are  described  in  Table  4. 


N 

s 

R 

C 

L 

N 

— 

N->S 

— 

N->C 

— 

S 

N->S 

S->R 

S->R 

S->R 

s->c 

c->s 

— 

R 

— 

S->R 

— 

— 

— 

C 

N- 

>C 

S->C 

C->S 

— 

C->L 

C->L 

C->L 

L 

— 

— 

— 

C->L 

— 

Table  4:  NSRCL  Interaction 


The  notation  used  in  this  Table  is  identical  to  that  used  in  Table  3  above.  Note  the 
interaction  when  an  active  supporter  meets  an  active  contrarian.  The  outcome  of  this 
meeting  (C  becomes  S,  S  becomes  C  or  no  change)  is  determined  by  the  probabilities  of 
conversion,  which  can  be  thought  of  as  the  relative  effectiveness  of  each  group’s 
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propaganda.  Note  also  that  when  either  the  parameters  or  initial  eonditions  are  sueh  that 
the  group  C  does  not  eome  into  play,  this  model  eollapses  to  that  proposed  by  Daley  and 
Kendall  (1964). 

The  system  of  differential  equations  for  this  model  follow: 


=-P^ns-p,nc 

at 

—=P,NS-P,S^-  P,SR  +  p,  SC 
dt 

P-=P,S-+P,SR 

dt 

—  =  P^NC  -P,SC-P,C^-P,CL 
dt 

^=P,&+li,CL 

dt 

with  initial  eonditions: 
iVo  =  P(l-(r,+r,)) 

So=CP 
R,  =  0 
Co=rcP 
D,=0 
r,+rc<l 


(5) 


The  definitions  for  the  parameters  in  this  model  are  similar  to  those  for  equations 
(1).  Note  that  all  yS'fsare  nonnegative  exeept  fory^^,  whieh  is  positive  if  the  aetive 
supporters  are  more  effeetive  than  the  aetive  eontrarians,  and  negative  otherwise. 


C.  CLOSED  FORM  SOLUTIONS:  THE  NSC  RECIRCULATION  MODEL 

In  this  seetion  we  present  two  models  whieh  are  speeial  eases  of  (2)  when  there 
are  no  latent  elasses;  both  aetive  supporters  (S)  and  aetive  eontrarians  (C)  may  beeome 
neutral  (TV).  We  eall  these  models  recirculation  models  beeause  there  are  no  absorbing 
states;  an  individual  may  ehange  elasses  infinitely  many  times.  These  models,  labeled 
NSC,  eonsider  the  flow  of  publie  opinion  as  opposed  to  ideology. 
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Both  NSC  models  asymptotically  converge  to  closed  form  solutions.  For  a 
detailed  description  of  the  methods  used  in  this  section,  see  Chapter  9  of  Boyce  and 
DiPrima  (1997). 

In  this  section  we  are  adjusting  our  notation  slightly.  Here  p^?,  represent 
recruitment  coefficients  into  the  S  and  C  classes  from  N,  while  s  represent  attrition 
coefficients  -  rates  of  flow  of  S  and  C  back  to  TV.  Additionally,  we  fix  P  =  N  +  S  +  C  =  \. 


1.  Quadratic  Model 

As  with  the  NSRC  and  NSRCD  models  we  assume  that  the  attrition  of  the  active 
spreaders  is  quadratic.  The  system  of  ODEs  representing  this  model  is; 

‘‘^^-p.NS-p.NC+pp+pp 

at 

P=p,NS-pp  (6) 

at 

^  =  p,NC-p,C'- 

This  model  has  the  interpretation  that  persons  become  politically  active  for  a  time, 
and  then,  after  their  period  of  activity,  they  simply  stop  caring  and  meld  back  into  the 
neutral  population  again.  Note  that  the  decay  term  for  both  the  supporters  and  the 
contrarians  is  quadratic.  The  quadratic  term  represents  the  following  phenomenon;  In  a 
large  population,  the  more  individuals  are  politically  active  on  a  given  side,  the  less  each 
individual  feels  he/she  must  contribute.  In  other  words,  the  rate  at  which  an  active 
opinionated  individual  (supporter  or  contrarian)  becomes  neutral  depends  on  how  many 
people  are  active  on  his  side;  the  more  active  people,  the  more  an  individual  would  be 
inclined  to  cease  his  activity. 

The  system  of  ODEs  in  (3)  has  a  single  asymptotic  stable  point  which  is  given  as 
a  function  of  the  model  parameters,  for  non- trivial  initial  conditions.  By  asymptotic 
stable  point,  we  mean  that  as  time  approaches  infinity,  the  size  of  each  class  will 
converge  to  a  single  point. 
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The  stable  point  of  this  system  ean  be  found  by  setting  the  derivatives  equal  to 
zero  and  solving  algebraieally.  Provided  that  none  of  the  /3^ '  s  are  equal  to  zero  and  the 
initial  eonditions  are  non-trivial  (there  are  initially  aetive  supporters  and  aetive 
eontrarians),  this  system  has  a  stable  point  at: 


N^  = 


S^  = 


C  = 


Pi 

A 

A2 

Pi 

A 

1  +  AL  + 

Pi 

A 

A 

Pi 

.  A 

1  +  AL  + 

P2 

A 

A 

(7) 


A  traee  of  the  system  for  given  parameters  and  initial  eonditions 
Ag  =.09,5q  =.01,Co  =.9  is  shown  in  Figure  6.  All  plots  in  this  Chapter  are  made  with 
MATLAB.  See  Seetion  D  for  details. 
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Figure  4.  Plots  of  the  NSC  Quadratie  model 


The  numerieal  eomputations  agree  with  our  theoretieal  asymptotie  results. 
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2,  Linear  Model 

In  the  second  model  we  assume  a  constant  attrition  rate;  fixed  proportions  of 


supporters  and  contrarians  become  neutrals.  The  system  of  ODEs  representing  this  model 
is: 


— =-p^NS-p^NC+p,S  +  p^  C 
dt 

^  =  p^NS-(i,S  (8) 

dt 

^  =  p,NC-fi,C 
dt 

The  linear  nature  of  this  model  makes  each  individual’s  political  activity  period 
independent  of  the  number  of  active  individuals  of  his  side. 

This  model  requires  a  different  approach  than  the  quadratic  one.  We  will  assume 
non-trivial  initial  conditions  and  proceed  by  taking  cases  on  the  parameter  values. 

Case  1:  fd.  >  p^,i  =  1,2. 

In  this  case  the  attrition  coefficients  are,  respectively,  greater  than  the  recruitment 
coefficients.  In  layman’s  terms,  there  is  a  net  outflow.  Next,  we  consider  the  equation 
corresponding  to  the  supporters.  The  case  of  the  contrarians  (third  equation)  is  treated 
similarly.  At  stability  we  have, 

^=S(A'p,-/?,)=0  (9) 

dt 


Because,  by  definition,  0  <  5  <  1 ,  the  sign  of  the  derivative  is  determined  by  the 
quantity  {Np^-/3^).  Because  0<A<1,  this  quantity  is  always  negative  in  this  case 

dS 

because  /3^>  p^.  Therefore,  for  any 5*  >0,  — <0  and  therefore  the  S  population  is 

dt 

decreasing  to  zero.  An  identical  argument  holds  for  the  C  equation. 

From  this  discussion  we  conclude  that  the  trajectory  of  the  system  is  always 
pointed  towards  the  point 

^oo=0,Q=0,N„=l  (10) 
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and  that  this  is  the  stable  point. 


Figure  5  depicts  the  trace  of  the  system  for  Case  1 . 
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Figure  5.  NSC  linear  model  for  Case  1 

In  Figure  5,  =  .5,4  =.3,Co  =.2,;?;  =.2,y0j  =.3,;72  =.5,y02  =-8 

Case  2:  >  p^y,p^  <  p^  <  p^ 

We  will  mirror  our  analysis  in  Case  1  and  examine  the  case 
where  P\>  P\,  P2  Pi  ’  realizing  that  the  result  is  symmetrical  for  the  other  case. 

Employing  the  same  reasoning  as  in  Case  1,  it  follows  that  the  S  population 
decreases  to  zero.  Therefore,  the  resulting  stable  point  is  only  a  function  of  the  N  and  C 
populations,  and  comes  from  solving  the  system  of  equations; 

CSP,-N^pP  =  CpN^p,-p,)  =  Q  (11) 

If  Coo  =  0 ,  then  the  last  equation  in  (5),  and  the  fact  that  Cg  >  0 ,  imply  that  for  t 
large  enough  (P2N  -  P2)  <  t) .  But  as  t  gets  larger  N  Since  P2  ^  Pi  h  follows  that 
{P2N- P2)  ^  0  ,  in  contradiction.  Therefore  ^  Oand  we  can  solve  (9)  algebraically  to 
obtain; 


^  _^_P2 


Pi  Pi 


(12) 
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By  symmetry,  the  case  where  ^2^  Pi')  P\  the  solution; 


'V„=^.c„=o,s„=i-^ 

A  Pi 


We  present  a  trace  of  this  system  as  Figure  6. 


(13) 


15 

time 


Figure  6.  NSC  linear  model  for  Case  2 
In  Figure  6,  =  .5,5'o  =.3,Co  =.2,;7i  =.2,y0j  =.3,;?2  =-5,y^2  =-2 

Case  3:  y0.  <  p^,  /  =  1,2 


This  case  has  two  sub  cases; 


(i) 


A 


Pi 

Pi 

And 

(ii) 

A 

Pi 

Pi 

Case  3i 

_A 


A 


It  is  immediate  that 5*^  x  =0  ,  since  otherwise  =—  and  N^=—  ,  which  is 

Pi  Pi 


a  contradiction.  Suppose  that; 


A<A<i.  (14) 

Pi  Pi 
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If  5'oo  =  Coo  =0  then  =  1 ,  and  therefore,  from  the  first  equation  in  (5),  it  follows 
that  for  t  large  enough  S{/3^  -Np^)  +  C{/32  -Np^)  >  0  must  hold.  But  for  t  large  enough 
N  ^  \  and  therefore  {/3.  -Np^)  <  0,/  =  1,2 ,  in  contradiction.  Hence  we  conclude  that  one 
of  the  following  solutions  is  correct: 

(a)iV„=A,c,=l-A,S.=0 

Pi  Pi 


Pi  Pi 

Suppose  (a)  is  true.  Let  £■  >  0  be  such  that 


s  < 


Pi 


Pi 


(15) 


The  right-hand-side  of  (13)  is  positive  because 

Pi  Pi 


Let  f  be  large  enough  such  that 


N 


(16) 


For  t>t'.  Such  a  t '  exists  because  N 


Pi 


From  (13)  and  (14)  it  follows  that 

Np^  -  >  0  for  >  L . 


Let  t”>t'  such  that 


A(Z")>0 


and 


A'(t")<0 


Such  Z"  exists  because  Sq>0  and  5*  ^  0. 


(17) 


(18) 


(19) 
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From  (15)  and  (16)  we  obtain  that  5"(t")  >  0 ,  which  contradicts  (17). 


Therefore  (b)  is  the  correct  solution  for  case  3i. 

By  symmetry,  for  the  case  where  —  <  —  <  1  the  solution  is  given  by; 

Pi  Pi 


h 


2_ 

Pi 


N  =—  /  =0  C  =1- 


Pi 


(20) 


We  present  a  trace  of  Case  3i  as  Figure  7: 
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Figure  7.  NSC  Linear  model  for  Case  3(i) 

In  Figure  7,  =  .5,5o  =.3,Co  =.2,y0j  =.3,;?2  =-5,  A 


Case  3ii 

From  the  second  and  third  equations  of  (5)  we  have  that 


N  = 


Spi  Pi 


c 

Cp, 


(21) 


Because 


A 

Pi 


A 

Pi 


it  follows  that  (19)  is  equivalent  to  the  differential  equation; 


p,[logSY  =  p,[logC]'. 


(22) 
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Since  and  cannot  be  both  zero,  it  follows  that 


N  = 


Pi 


K 

Pi 


and  from  (19)  and  (20)  we  obtain  as  a  solution  of; 


(23) 


f 

1 

V 


Pi 


A 

Pi 


c  Pi 
P\ 

^0 


(24) 


Since  =1 ,  we  obtain  also  a  value  for  .  We  present  a  trace  of  case 

3(ii)  as  Figure  10. 
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Figure  8.  NSC  linear  model  for  Case  3(ii) 

In  Figure  8,  =.5,5*0  =.3,Co  =.2,;?;  =.4,y0j  =.3,;72  =-8,  A  =-6 

D.  NUMERICAL  ANALYSIS 
1,  Numerical  Methods 

In  Section  C  we  developed  closed-form  solutions  to  some  special  cases  of  the 
general  models.  Analytic  solutions  to  the  cases  that  we  are  most  interested  in,  namely 
NSRL  and  NSRCL,  are  more  elusive.  We  proceed  with  numerical  analysis  of  these  two 
models. 


24 


The  fact  that  we  do  not  possess  a  closed-form  solution  to  the  cases  of  interest  to 
us  is  not  as  great  a  difficulty  as  it  may  initially  seem.  We  may  approximate  the  solution 
to  the  ODE  system  numerically.  For  solving  the  two  aforementioned  models  and  perform 
the  analysis  we  use  the  routine  ode45  supplied  with  the  software  package  MATLAB. 
The  basis  of  this  routine  is  the  4*  order  Runge-Kutta  method.  Interested  readers  are 
directed  to  Chapter  8  of  (Boyce  and  DiPrima,  1992). 

For  the  NSRL  model  we  assume  a  base  case  where  all  the  transmission 
coefficients  are  equal  to  1 .  As  before,  we  assume  that  P=\. 


2.  The  NSRL  Model 

For  the  NSRL  model  we  assume  a  base  case  where  all  the  transmission 
coefficients  are  equal  to  1 .  As  before,  we  assume  that  P  =\.  We  first  present  in  Figure 

9  the  base  case  trace  in  time  with  =  .  1 . 
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Figure  9.  NSRF  model;  base  case 


Notice  that  after  10  time  units  the  active  supporters  are  about  10%  of  their  initial 
number.  The  proportions  in  the  population  of  latent  supporters  (if)  and  latent  contrarians 


(L)  converge  in  the  long  run  to  .42  and  .34,  respectively. 
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We  next  investigate  the  effect  of  the  initial  sizes  of  the  S  (active  supporters)  andZ 
classes  on  the  composition  of  the  population  at  stability.  Specifically,  for  a  given  initial 
proportion  rs  of  active  supporters  we  ask;  what  is  the  initial  proportion  of  latent 
contrarians  ri  such  that  at  stability  the  number  of  latent  supporters  is  equal  to  the 

number  of  latent  contrarians  ?  This  Measure  of  Effectiveness  (MOE)  indicates  the 

“neutralization”  effect  of  the  two  classes  of  population  -  the  initial  conditions  that  lead  to 
parity.  Many  other  MOEs,  including  time-dependant  ones,  could  be  developed  and 
implemented. 

Clearly,  the  final  ratio  between  the  supporters  and  the  contrarians  is  a 

monotone  decreasing  function  of  the  initial  fraction  Tc  of  contrarians. 

We  use  a  binary  search  in  MATEAB  to  iteratively  approach  the  critical  proportion 
of  initial  contrarians  for  each  initial  proportion  of  pro-spreaders.  This  is  a  technique  we 
will  use  frequently  in  the  following  analyses.  We  set  the  effectiveness  parameters  at  their 
base  case  values.  Eigure  12  presents  the  proportion  of  initial  latent  contrarians  ri  required 
for  achieving  parity  for  a  given  proportion  rs  of  initial  supporters. 
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Eigure  10.  Proportion  of  initial  latent  contrarians  ri  required  for  achieving  parity  for  a 

given  proportion  rs  of  initial  supporters. 
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Figure  10  is  interesting  because  it  shows  that  for  the  given  modeling  assumptions, 
a  somewhat  large  increase  in  the  S  population  is  offset  by  a  modest  increase  in  the  L 
population. 

Recall  that  in  Figure  9,  the  N  population  (neutrals)  does  not  go  to  zero  as  time 
goes  to  infinity,  but  rather  reaches  a  stable  point.  We  next  consider  the  size  of  the  TV  class 
at  parity  as  a  function  of  the  initial  proportion  of  supporters  r^.  Note  that  the  parity 
condition  uniquely  determines  the  value  of  ri,  as  shown  in  Figure  10. 
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Figure  1 1 .  N(infinity)  at  parity  for  increasing  numbers  of  initial  supporters 

Figure  11  shows  that  the  number  of  neutrals  approaches  zero  rather  rapidly  at 
parity,  as  the  initial  number  of  supporters  increases.  A  surprising  result  from  this  graph  is 
that  the  introduction  of  class  L  causes  A( infinity)  to  be  larger  than  .203,  which  is  the 
expected  proportion  from  Daley  and  Kendall  (1964).  We  attribute  this  to  the  fact  that  the 
attrition  of  class  S  is  faster  in  our  model  than  in  theirs,  making  fewer  available  to  recruit 
out  of  class  N. 

We  next  consider  the  implications  of  varying  the  effectiveness  parameter  of  the 
latent  contrarians,  .  We  vary  the  parameter  and  find  the  value  of  the  parameter 

that  is  required  to  achieve  parity.  We  fix  =.l  throughout.  Figure  14  shows  the  result 
of  this  experiment. 
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Figure  12.  Critical  with  varying  effectiveness (3^  ,rs=.\ 

Note  in  Figure  12  that  as  ^04  approaches  1,  approaches  .102,  which  is  the  value 
shown  in  Figure  12  for  =  .  1 . 

We  next  consider  the  behavior  of  the  model  when  we  vary  the  parameters  (3^  and 
/3^ ,  which  are  the  S  recruitment  rate  from  the  N  population  and  the  L  conversion  rate 
from  the  S  populations,  respectively.  We  fix  the  values  of  and  at .  1  each  (as  in  the 
base  case)  and  search  for  parameter  values  that  achieve  parity. 
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Figure  1 3 .  Contrarian  effectiveness  required  to  achieve  parity,  =  /y  =.  1 . 
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Note  that  as  the  conversion  effectiveness  increases,  a  very  small  increase  in  the 
effectiveness  of  the  contrarians  must  be  countered  by  a  very  large  increase  in  the 
effectiveness  of  the  active  supporters  in  order  to  attain  parity.  For  example,  when 
^^4  =  1,  then  =  0.75  at  parity.  Flowever,  when  ^^4  =  2,  then  =  7.5  at  parity,  and 

y^4  =  2.2,  results  in  =  15 . 

Our  numerical  analysis  has  focused  on  the  scaled  version  of  the  process.  We  now 
consider  the  implications  of  using  an  un-scaled  process,  allowing  the  population  size  P  to 
grow  while  leaving  the  number  of  initial  active  supporters  So  fixed.  This  mirrors  the 
analysis  performed  by  Daley  and  Kendall  (1964).  In  Figure  16,  we  fix  the  number  of 
initial  spreaders  at  1,  and  allow  the  population  size  to  grow. 
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Figure  14.  Critical  ,  5*0=1 

As  the  population  increases,  the  proportion  of  contrarians  required  to  achieve 
parity  appears  to  be  invariant  to  the  population  size.  The  source  of  the  variability  in  Fig. 
14  is  an  artifact  of  double-precision  arithmetic  and  not  a  property  of  the  system  itself. 

The  deterministic  model  has  the  behavior  that  the  final  distribution  of  classes  N, 
R,  and  L  is  invariant  to  a  fixed  5^  for  large  populations.  Stochastic  models  do  not  exhibit 

this  behavior  and  is  an  example  of  a  case  where  deterministic  models  are  unsatisfactory 

for  analysis.  We  address  this  issue  in  the  following  chapter. 
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3.  NSRCL 

Next  we  eonsider  the  symmetrie  model,  whieh  has  more  variables  than  the  NSRL 
model.  For  the  NSRCL  model  we  assume  a  base  case  where  all  the  transmission 
coefficients  yS'.  are  equal  to  1  except  13^=0  (recall  that  is  governs  the  ‘conversion’  of 

active  members  of  the  pro-  and  counter-  classes,  an  asymmetric  behavior).  As  before,  we 
assume  that  P  =\. 

Figure  17  below  presents  the  trace  of  the  model  for  the  base  case. 


NSRCL  Base  Case 


Figure  15.  NSRCL  Base  Case 


Note  that  due  to  the  symmetric  nature  of  the  model,  the  traces  of  S  and  C  as  well 
as  R  and  L  are  coincident. 

We  next  present  four  time  traces  of  the  NSCRL  model,  slightly  perturbed  from 
the  Base  Case,  in  order  to  show  the  behavior  of  the  population  classes  (Class  N  is  not 
shown  for  the  sake  of  clarity).  We  show  these  plots  as  a  contrast  to  Figure  15,  showing 
the  behavior  of  the  classes  under  non-homogeneous  parameter  values. 
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rS  =  .1,  rC=  .1,  p1  =  .3,  p2  =  .5 


rS  =  .1,  rC=  .15,  pi  =  .3,  p2  =  .5 
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rS  =  .2,  rC  =  .2,  pi  =  .3,  p2  =  .5 
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Figure  16.  NSRCL  Model,  with  varying  initial  conditions 

Similar  to  our  analysis  of  the  NSRL  model,  we  next  consider  the  proportion  of 
spreaders  of  each  type  to  achieve  parity  in  the  base  case  (all  conversion  intensities 
equal  to  1). 
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Figure  17.  Proportion  of  active  contrarians  to  achieve  parity  for  a  given  initial 

proportion  of  active  supporters. 
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As  the  NSRL  model  is  symmetric,  the  linear  plot  in  Figure  17  is  not  surprising. 
We  repeat  the  analysis  for  a  slightly  different  case.  Here  we  fix  the  parameter  at  .5, 
making  the  contrarians  half  as  effective  at  recruitment  from  the  neutral  population  as  the 
supporters. 
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Figure  18.  Proportion  of  active  contrarians  Vc  to  achieve  parity  for  a  given  initial 

proportion  of  active  supporters,  =  .5 

From  Figure  18  we  see  that  the  effect  of  lower  contrarian  recruitment 
effectiveness  depends  on  the  proportion  of  supporters  ;  for  small  values  of  this 
parameter  the  critical  proportion  of  initial  contrarians  is  considerably  higher  than  , 
while  for  larger  values  of  the  behavior  is  similar  to  the  case  of  ^  shown  in  Figure 
19. 

We  conclude  our  analysis  with  respect  to  y^j^Y  fixing  the  initial  proportion  of 
supporters  at  its  base  case  {r^  .1)  and  varying  the  effectiveness  of  recruitment  by  the 
contrarians  y^j .  Figure  19  presents  the  results  of  this  analysis. 
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Figure  19.  Critical  proportion  of  contrarians  to  achieve  parity  for  varying  contrarian 

recruitment  effectiveness  [3^ .  r^.\ 

The  results  in  Figures  18  and  19  are  surprising  at  first;  we  might  expect  parity  to 
be  unachievable  for  low  recruitment  effectiveness.  The  explanation  for  this  phenomenon 
is  that  the  ability  to  recruit  is  offset  by  having  such  a  large  proportion  of  the  population 
on  the  counter  side  initially. 

We  repeat  the  analysis  in  Figure  19,  but  instead  of  varying  ,  we  set  =  /3.j  and 

vary  these  parameters  together.  Recall  that  these  parameters  govern  the  rate  of  conversion 
for  the  contrarians,  C,  into  the  absorbing  state  L  both  via  C-C  interaction  and  C-L 
interaction. 
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Figure  20.  Critical  proportion  of  Class  C  to  achieve  parity,  varying  =  P-j,  all  other 

parameters  as  in  the  base  case 

An  interpretation  of  varying  is  that  for  small  values,  the  active  life  of  an 

active  Contrarian  is  very  long;  for  large  values,  it  is  very  short.  We  see  here  that  as  their 
life  span  increases  towards  infinity,  the  proportion  required  to  achieve  parity  becomes 
very  small.  As  their  lifespan  becomes  very  short  (large  values  of  the 

proportion  required  approaches  .4,  which  is  to  say  that  their  recruitment  is  virtually  nil. 

The  NSRCL  model  has  no  invariance  property  analogous  to  the  behavior  shown 
in  Figure  14. 

The  final  issue  we  consider  is  the  effect  of  varying  the  number  of  initial 
supporters,  initial  contrarians,  and  the  effectiveness  of  conversion.  Recall  that  is  the 

rate  at  which  supporters  ‘convert’  contrarians.  This  is  the  only  parameter  in  our  model 
that  is  allowed  to  take  on  negative  values  (a  negative  value  means  that  contrarians 
convert  supporters). 

We  vary  the  parameter  along  the  X  axis  of  Figure  22;  we  vary  along  the  Y 
axis.  The  z  axis  is  .  We  fix  = .  1  throughout. 
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Because  the  model  is  symmetric,  the  same  results  would  be  obtained  if  everything 
were  reversed;  that  is,  what  we  consider  to  be  supporters  and  contrarians  in  this  analysis 
is  unimportant  as  long  as  we  are  consistent.  In  other  words,  we  could  have  just  as  well 
varied  and  /3^ ;  the  plot  obtained  would  be  a  mirror  image  of  the  ones  shown  below. 


Conversion  effectiveness 


Figure  21.  The  numerical  difference  between  and  as  a  function  of  model 

parameters  (four  views) 

Figure  21  presents  four  views  of  the  same  plot.  We  see  from  looking  at  the  plot 
oriented  along  the  axis  that  in  the  ranges  considered  here,  the  initial  conditions  have  a 
greater  effect  than  the  effectiveness  of  conversion. 
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IV  STOCHASTIC  MODELS 


A.  MOTIVATION 

While  in  Chapter  III  we  presented  a  deterministie  approaeh  to  the  ideologieal 
propagation  problem,  in  this  ehapter  we  obtain  additional  insights  by  studying  the 
stoehastie  proeess  that  deseribes  the  ideologieal  propagation  among  different  population 
eategories.  The  main  benefit  of  stoehastie  models  as  opposed  to  deterministie  ones  is  that 
they  enable  the  analysis  of  random  interaetions  among  individuals;  the  eost  is  that  the 
resulting  models  beeome  less  traetable. 

Beeause  NSC  models  laek  absorbing  states,  the  teehniques  employed  in  this 
ehapter  are  not  appropriate  for  the  NSC  reeireulation  models  introdueed  in  Chapter  III.  In 
this  ehapter  we  foeus  on  the  asymmetrieal  model  (NSRL).  Reeall  that  the  NSRL  model 
has  four  different  types  of  aetors: 

•  N:  the  neutrals,  who  are  members  of  the  population  who  earry  neither  the  pro-  or 
eounter-  message. 

•  S:  the  aetive  supporters,  who  earry  the  pro-  message  and  are  aetively  spreading  it. 

•  R:  the  retired  supporters,  who  earry  the  pro-  message  but  are  not  spreading  it. 

•  L:  the  latent  eontrarians,  who  earry  the  eounter-  message  but  are  not  aetively 
spreading  it. 

The  engagement  rules  among  these  groups  are  the  same  as  those  proposed  in 
Table  3  of  Chapter  III.  One  differenee  is  that  the  parameter  rj  in  our  stoehastie 

development  takes  the  plaee  of  the  parameter  in  the  deterministie  model;  rj  deseribes 

the  effeetiveness  of  the  Latent  Contrarians  L  in  eonverting  members  of  elass  S.  We 
ehoose  this  as  the  single  degree  of  freedom  in  the  parameters. 

The  engagement  rules  deseribed  in  Table  3  imply  that  the  proeess  is  strietly 
evolutionary,  so  that  an  absorbing  state  is  reaehed  in  finite  time  with  probability  one.  This 
faet  indieates  that  eapturing  the  main  properties  of  the  distribution  of  the  state  of  the 
proeess  when  it  enters  the  absorbing  states  -  the  states  from  whieh  no  further  dynamies 
oeeur  -  is  erueial  towards  understanding  the  implieations  of  our  models.  This  is  the  foeus 
in  this  ehapter. 
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This  chapter  is  organized  as  follows:  In  seetion  B  we  formulate  our  model. 
Seetion  C  presents  our  results  for  the  small  population  seenario,  while  Seetion  D  does  the 
same  but  for  large  populations.  Seetion  E  presents  the  eonelusions  of  this  Chapter. 


B,  MODEL  FORMULATION 

In  this  seetion  we  formulate  the  eontinuous  time  Markov  ehain  and  the  assoeiated 
diserete  time  ehain  that  form  the  eore  of  our  model. 

A  population  of  size  P  is  split  among  4  mutually  exhaustive  groups:  Neutrals  (N), 
Supporters  (S),  Retired  Supporters-  (R)  and  Latent  Contrarians  (L).  We  assume  that  the 
time  dynamies  of  these  groups  follow  a  eontinuous  time  Markov  ehain 

((A,5,R,L)(0:t>0),  (1) 

with  state  spaee 

S  =  {{n,s,r,  1)  e  :n  +  s  +  r  +  l  =  P,n>0,s>0,r>0,l  >0) ,  (2) 

and  initial  eonditions 

{N,S,R,L){0)  =  {P-S,-L„S„0,L,),  (3) 


where  and  Lq  are  non-negative  integers.  In  our  notation,  N(t)  is  the  number  in  elass  N 
at  time  t;  S{t),  R{t)  and  L{t)  are  likewise. 

Given  that  for  all  t>0  we  have  R(t)  =  P - N(t)- S{t)- L(t)  ,  mueh  of  what 
follows  is  without  loss  of  generality  framed  in  terms  of  the  (A,5',L)(.)proeess,  whieh  for 
notational  eonvenienee  we  define  as  'K  =  ((N,S,L)(t):t>0) .  The  transition  rates  are 
given  by 


(n, s,l)^(n-\,s  + 1, /)  at  rate  13ns 
for  the  N-S  interaetion,  and  by 


{n,  s,  /)  -^{n,s-  2,  /)  at  rate  /3 


(4) 

(5) 


for  the  S-S  interaetion.  For  the  S-R  interaetion  the  transition  rate  is 
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{n, s,  1)  -^{n,s- 1, /)  at  rate  Ps{P -n-s-l). 

The  S-L  interaetions  are  governed  by  the  transition  rate 

{n,s,l)  ^  {n,s  -  \,l  +  \)  at  rate  Prjsl.  (6) 

The  parameter is  the  stoehastie  mixing  rate  parameter  and  ri,0<ri<\,  is  the 
effeetiveness  of  the  eontrarians.  Figure  22  illustrates  the  possible  transitions  out  of 
state  (u,  5,/) . 


Figure  22.  State  Transitions 


The  evolutionary  nature  of  the  proeess  X ,  as  evideneed  by  its  transition  rates, 
implies  that  Xreaehes  an  absorbing  state  in  finite  time  with  probability  1.  Flenee,  in 
order  to  study  the  distribution  of  the  absorbing  state  it  sufliees  to  eonsider  the  assoeiated 
diserete  time  Markov  ehain  (DTMC)  ofX  .  For  this  reason  let 

Y.=((Af,5,i),:»>0)  (7) 

be  a  diserete  time  Markov  ehain  defined  on  the  same  state  spaee  and  with  the  same  initial 
eonditions  as  X  .  Due  to  finite  population,  for  all  u  >  0  the  relationship 

R  =P-N-S-L  (8) 

n  n  n  n  v  / 
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holds,  so  we  know  the  amount  of  the  retired  pro  spreaders  by  storing  information  about 
only  three  groups.  We  now  consolidate  equations  (1-8)  and  derive  the  transition  matrix 
for  Y: 


{N,S,L)„^,={N,S,LX  + 


(-1,+1,0)  w.p.  (a) 

J  nsl 

(0,-2,0)  w.p.  {b) 

J  nsl 

(0,-l,0)  w.p.  ^  ^  ^  ^  (c) 

J  nsl 

(0,-1, +1)  w.p.  ^  {d) 

J  nsl 

Where 

L,=p-y2s+vi-i-y2  (10) 

The  transition  matrix  is  conditional  on  a  transition  occurring;  we  scale  the  rates  of 
transition  in  equations  (9)  by  the  factor  in  equation  (10)  to  have  a  matrix  where  the 
probabilities  sum  to  1 .  Additionally,  the  rate  parameter  /3  has  divides  out. 

These  transition  rules  indicate  that  the  absorbing  states  of  Y  are  given  by 

B  =  iin,s,l)eS:s  =  0),  (11) 

and  that  P{t  <  oo)  =  1 ,  where 


r  =  inf{n>0:Y^  gB}  (12) 

Due  to  the  evolutionary  nature  of  the  process,  we  can  place  a  firm  deterministic 
bound  on  the  value  of  r  .  Consider  the  limiting  case  on  the  initial  conditions  where  the 
population  consists  of  1  Y  and  P-\  N.  Each  member  of  the  susceptible  population  must 
pass  through  class  S  en-route  to  an  absorbing  state.  Now  consider  a  (possible  but 
unlikely)  limiting  case  where  N^=0.  Each  member  of  the  susceptible  population  must 

pass  through  class  S  en-route  to  an  absorbing  state,  for  E’-1«E’ transitions.  The 

maximum  number  of  transitions  out  of  class  S  occurs  when  S-S  stifling  does  not  occur. 
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This  path  is  unlikely,  but  possible.  If  we  eonsider  the  maximum  number  of  transitions 
out  of  class  S,  we  find  it  is  P.  Therefore,  r  <  2P .  Note  that  this  is  a  limiting  case  in  the 
absence  of  contrarians.  The  addition  of  contrarians  will  necessarily  reduce  the  maximum 
number  of  transitions  further. 

C.  RESULTS  FOR  SMALL  POPULATIONS 

1.  Algorithm 

Equations  (9)  and  (10)  contain  all  the  information  about  the  process  Y,  but  some 
care  is  required  to  extract  useful  information  about  the  distribution  of  of  the  state  of  the 
process  when  it  enters  the  absorbing  states;  that  is  the  distribution  of  the  number  of 
individuals  in  classes  R  and  L  when  the  number  of  individuals  in  class  S  becomes  0.  A 
naive  direct  approach  would  involve  mapping  the  three-dimensional  state  space  onto  one 
dimension,  but  it  quickly  becomes  intractable  for  non-trivial  population  sizes,  and  the 
resulting  transition  matrix  is  very  sparse,  which  leads  to  numerical  issues  due  to 
insufficient  computer  memory. 

A  more  refined  approach  that  is  numerically  stable  for  population  sizes  in  the 
order  of  the  hundreds  is  to  proceed  sequentially,  as  described  next. 
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Algorithm:  Compute  Absorbing  Distribution 
Initialize : 


PR. 


-{p+i)(i=+i)(p+i) 


TEMP 


(i=+i)(p+i)(p+i) 


Absorbing^p^.^^p^,^ 


=  0 
=  0 
=  0 


while  i^any  element  of  PR>0){ 


TEMP  =  PR 


choose  one  element  of  PR  to  updatej 

ApplyEqn(9)to  Pi?  to  determine  successor  state  probabilities. 


Add  to  TEMP 


successor  State 


if  {successor  stateis  Absorbing)  { 
Add  to  Absorbing 


Set  (TEMf 


successor  State 


=  0) 


} 

PR  =  TEMP 


} 


Algorithm  1 


The  3-dimentional  matrix  PR  contains  the  probabilities  associated  with  each  non- 
absorbing  state.  It  has  dimension  P  +  \  in  order  to  include  zero  as  a  value  for  each 
population. 

Each  State  {N,S,L)  has  up  to  four  successor  states  and  multiple  predecessor 
states;  a  ‘breadth-first’  approach,  such  as  above,  is  required.  It  is  useful  to  think  of  the 
analogy  of  a  wave  moving  out  from  a  stone  dropped  in  a  pond.  This  is  in  essence  how  the 
probability  mass  travels  from  the  (given)  initial  state.  The  non-zero  elements  of  PR  are 
on  the  ‘crest’  of  the  wave.  The  successor  probabilities  are  computed  step-wise;  when  the 
successors  of  any  given  state  are  computed,  the  wave  has  passed.  Algorithmically,  that 
state  is  no  longer  considered;  its  probability  is  assigned  to  zero.  This  process  continues 
until  all  of  the  probability  mass  has  diffused  to  absorbing  states. 

TEMP  is  merely  a  temporary  copy  oIPR. 
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Absorbing  is  a  (P  +  l)x(P  +  l)  matrix,  holding  the  probability  mass  associated 
with  the  number  of  individuals  in  classes  R  and  L  when  the  number  of  individuals  in 
class  S  becomes  0,  defined  in  equation  (11).  Note  that  because  S  is  zero  for  the 
absorbing  states,  two  dimensions  are  sufficient  to  contain  the  distribution  of  X. 

At  each  jump  epoch,  the  probability  mass  diffuses  from  each  state,  eventually 
reaching  an  absorbing  state  where  it  accumulates.  When  the  process  is  ‘near’  a  boundary 
(i.e.  one  jump  from  a  disallowed  state),  some  care  must  be  exercised  to  ensure  that  the 
process  remains  in  the  non-negative  integers.  The  only  boundary  rule  that  needs  to  be 
explicitly  enforced  is  to  disallow  the  S  to  S  stifling  process  when  only  one  spreader  is 
active.  We  implement  the  above  in  MATLAB,  with  the  code  found  in  Appendix  1. 

2,  Numerical  Results  and  Analysis 

We  report  the  results  of  our  computations  with  respect  to  the  distribution  of  the 
number  of  individuals  in  classes  R  and  L  when  the  number  of  individuals  in  class  S 
becomes  0.  The  only  limit  on  the  accuracy  of  the  computation  is  the  numerical  precision 
of  the  computer.  We  show  the  plots  below  for  four  cases: 
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Figure  23. 


Case  I;  P=100,  4  =1,  =1, 77=1 


0.04  V- 


0.03- 


Number  in  class  R  0  0  Number  in  class  L 

Figure  24.  Case  II :  F’=I00,  4  =!>  5'o  =!>  ^  =  -5 


0.025 


Number  in  class  R  Number  in  class  L 


Figure  25.  Case  III; P =100,  4=10,  4=10,  ri  =  .\ 
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X  10 


Number  in  class  R  °  0  Number  in  class  L 

Figure  26.  Case  IV; /*= 100,  Lq=\0,  50=10,  ri=\ 

In  Figures  23  and  24  5o=l  and  there  are  positive  limiting  probabilities  for  the 

number  of  individuals  in  elass  R  and  class  L  being  close  to  0.  The  sum  of  these 
probabilities  is  the  probability  that  the  ideological  spread  ‘fizzled  out’  and  does  not  gain 
traction  throughout  the  population  as  a  whole;  that  is,  the  population  remains  neutral. 

Figures  25  and  26  show  cases  where  the  process  spreads  through  the  population. 
These  are  in  stark  contrast  to  Figures  23  and  24,  where  there  is  a  large  ‘spike’  near  the 
origin.  This  spike  represents  the  probability  that  nothing  significant  happens. 

The  probabilities  displayed  in  Figures  25  and  26  look  characteristically  ‘bell¬ 
shaped’  or  Gaussian;  in  the  next  section  we  argue  that  they  are  approximately  of  that 
form. 

D.  RESULTS  FOR  LARGE  POPULATIONS 

The  goal  in  this  section  is  to  study  the  distribution  of  the  number  of  individuals  in 
classes  R  and  L  when  the  number  of  individuals  in  class  S  becomes  0  for  large  population 
sizes,  a  situation  where  Algorithm  1  cannot  be  implemented  due  to  memory  constraints. 
As  we  shall  momentarily  see,  the  analysis  of  this  section  lends  support  to  the  apparent 
normality  of  the  distribution  of  the  number  of  individuals  in  classes  R  and  L  at  the  time 
the  number  of  individuals  in  class  S  becomes  0  observed  in  Figures  25  and  26.  The  basis 


45 


for  this  section  are  the  papers  Barbour  (1972)  and  Barbour  (1973),  who  proves  a 
functional  Central  Limit  Theorem  for  Markov  processes  of  the  type  described  in  this 
chapter,  and  provides  an  analytically  tractable  approach  to  solve  the  differential  equations 
that  arise  in  this  setting.  The  focus  on  this  section  is  on  applying  the  results  of  these  two 
papers  to  our  propagation  of  ideas  problem.  The  key  step  is  to  normalize  the  population 
size,  and  relax  the  requirement  for  an  integer  state  space.  In  this  section,  we  consider  the 
proportions  (non-integer)  of  the  population  in  each  class  R  and  L  as  opposed  to  the 
(integer)  members  in  each  class. 

1.  Theoretical  Results 

With  this  in  mind,  consider  a  sequence  of  CTMC’s  (X^)  indexed  by  the 
population  parameter  P,  and  defined  on  t  >  0 .  Each  process  Xp  has  initial  condition 

{N,S,L){0)  =  {P{\-r,-rJ,Pr„Pr,)  (13) 

for  some  constants  >  0,  r^>0,  r^+r^  <  \ .  Because  the  initial  size  of  each  class  grows 

with  the  population  size,  our  initial  conditions  suggest  that  the  propagation  will  not  die 
out  in  its  initial  stages,  as  observed  in  Figures  23  and  24. 

To  complete  the  framework,  let 

Xp(^)  =  E 'Xp(0  for  0  <  t  <  T  (14) 


be  a  sequence  of  stochastic  processes,  where  each  Xp(t)  represents  the  proportion  of 
people  in  a  population  of  size  P  that  at  time  t  belong  to  each  class  N,  S,  ox  L. 

The  transition  rates  for  the  scaled  process  Xp  are  given  by 
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subject  to  the  boundary  conditions 


n>0,  5>0,  /  >0  (16) 

It  can  be  verified  that  \p{t)  meets  the  conditions  needed  to  apply  the  results  of  Barbour 
(1972)  and  (1973).  The  result  in  these  papers  that  is  most  useful  for  our  purposes  is  that 

sup,  \P{Zp{t)  <  t;)-O(t;/E(0)|  =  0(P  ^'^  logP)  (17) 

for  fixed  t,  0<t<T ,  where: 

0(.)  is  the  distribution  function  of  a  standard  normal  variable; 

E(t)  is  a  deterministic  variable  that  depends  on  t  and  the  parameters  of  the  process; 

Zp{t)  =  P'^'^  {xp(t)-^(t)}*/2(^(t),t)  is  the  normalized  version  of  Xp(t),  where 

^{t)  is  the  deterministic  solution  found  in  Section  2  of  Chapter  III 

h(,.)  is  a  deterministic  vector-valued  function  that  depends  on  the  problem  parameters; 

T  is  the  largest  time,  possibly  infinite,  for  which  we  can  find  a  deterministic  solution 
utilizing  the  methods  in  Section  2  from  Chapter  III. 

The  result  is  that  a  normalized  version  of  Xp(t)  converges  in  functional  space  to  a 
Normal  random  variable  with  variance  that  depends  on  the  initial  conditions  and  the  rate 
parameters,  asP  ^  oo ;  in  words,  this  means  that  the  distribution  ofZp(t)is  never  further 

than  0{^P  logf^  from  the  Normal  distribution.  This  result  is  operationally  significant. 

First  it  implies  that  the  transient  regime  and  the  state  upon  absorption  of  the  stochastic 
model  that  describes  the  propagation  of  ideas  is  random.  Secondly,  for  any  time  t, 
0  <  t  <  r ,  the  distribution  of  the  proportion  of  individuals  in  each  ideological  group  is 
approximately  Normal  with  mean  given  by  the  deterministic  solution,  and  covariance 
structure  that  can  be  analytically  computed  and  depends  on  the  initial  conditions  and  rate 
parameters. 

In  the  next  subsection  we  illustrate  these  results. 
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2,  Numerical  Results  and  Analysis 

In  this  subsection  we  present  various  illustrative  scenarios.  We  start  by  comparing 
the  mean  of  stochastic  (found  using  Algorithm  1)  with  the  deterministic  solution. 
Specifically,  we  are  interested  in  the  quantity; 


X 


00,  deterministic 


-E{Y^ 


r, stochastic  -I 


P 


(18) 


where  t  is  as  defined  in  eqn(12),  for  large  values  of  P.  In  the  following  plots,  the 
parameter  rj  (stochastic)  is  set  to  1,  as  well  as  the  parameters  ^0'.  (deterministic).  Note 

that  although  we  only  have  theoretical  grounds  for  expecting  the  means  to  converge  when 
the  initial  conditions  are  posed  as  a  proportion  of  the  population  P,  this  behavior  is  also 
exhibited  for  fixed  numbers  as  the  population  grows.  We  expect  the  convergence  to  be 
monotone  in  the  cases  where  the  problem  is  posed  as  proportions;  Case  II  shows  non¬ 
montonicity  which  may  result  if  the  initial  conditions  are  posed  as  fixed  numbers. 

CASE  I;  is  set  at  1 ,  and  is  set  at  1 
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Figure  27. 
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CASE  II;  5o  is  set  to  1,  and  is  set  to  10%  of  the  population 
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Figure  28.  Case  11  eonvergence  of  Means 


CASE  III;  Sq  is  set  to  10%  and  is  set  to  1 
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Figure  29.  Case  III  convergence  of  Means 


CASE  IV;  is  set  to  1  and  Zg  is  set  to  30%  of  the  population 
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Figure  30.  Case  IV  convergence  of  Means 


In  each  case  presented  here,  the  quantity  defined  in  equation  (18)  tends  toward 
zero  with  increasing  population  size,  as  expected.  We  feel  confident  in  stating  that  our 
numerical  results  for  means  does  agree  with  Barbour  (1972)  and  (1973)  papers. 

We  now  illustrate,  through  various  choices  of  initial  conditions  and  rate 
parameters,  the  fact  that  the  distribution  of  the  number  of  individuals  in  classes  R  and  L 
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at  the  time  the  number  of  individuals  in  class  S  becomes  0  is  well  approximated  by  a 
multivariate  normal  distribution  as  the  population  size  increases. 

Now  we  consider  a  run  of  the  model  with  the  parameters  P  =  100,?7  =  1.  The 
bivariate  distribution  of  the  number  of  individuals  in  classes  R  and  L  when  the  number  of 
individuals  in  class  S  becomes  0  under  different  initial  conditions  are  plotted  here; 


Figure  31.  F’=100,5'o=10,To=10,;7  =  l 
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Figure  32.  P  =  100,S^f=20,L^^=20,r/  =  l 
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Figure  34.  F’=200,5'o  =15,Zo  =20,;;  =  1 

These  plots  illustrate  that  for  some  initial  conditions  as  specified  in  Equation  (13), 
the  distribution  of  the  number  of  individuals  in  each  class  R  and  L  when  the  number  of 
individuals  in  class  S  becomes  0  is  approximately  multivariate  Normal.  However,  note 
the  apparent  multimodal  behavior  displayed  in  figure  33.  We  make  this  assertion  based 
on  the  theoretic  framework  in  Barbour  (1972  and  1973).  In  essence,  we  are  applying  an 
extension  of  the  Central  Limit  Theorem  to  study  the  expected  value  of  a  more  general 
Markov  process. 

Barbour  (1972)  provides  a  procedure  for  computing  the  variance/covariance 
structure  of  the  limiting  Gaussian  distribution  of  proportion  of  individuals  in  each 
absorbing  class.  An  application  of  this  method  is  beyond  the  scope  of  this  thesis. 
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E.  SUMMARY 

Our  results  bring  home  that  for  the  areas  of  most  interest  for  persuading/ehanging 
publie  opinion,  the  deterministie  model  represents  the  mean  outcome  of  the  stochastic 
propagation  process.  In  the  situation  where  the  ideological  spreading  does  not  die  out  in 
its  initial  stage,  the  results  of  this  chapter  indicate  that  the  final  proportion  of  people  in 
each  of  the  ideological  classes  is  approximately  Normally  distributed  with  a  mean  given 
by  the  solution  of  a  system  of  ODEs  similar  to  those  developed  in  Chapter  III  and  a 
covariance  structure  that  can  be  computed  from  the  rate  parameters  and  the  initial 
conditions. 

We  believe  that  the  major  contribution  of  this  chapter  is  to  provide  a  framework 
to  explore  models  of  ideology  propagation.  The  interaction  rules  chosen  for  this 
particular  development  are  reasonable  and  fit  in  smoothly  with  other  models  discussed  in 
the  literature. 
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V.  CONCLUSION 


A.  RESULTS 

In  this  thesis  we  have  extended  the  models  of  rumor  propagation  to  the  case  of  a 
competitive  ‘hearts  and  minds’  campaign.  We  have  developed  a  general  methodology 
and  applied  it  to  study  two  cases. 

In  Chapter  III  we  develop  two  deterministic  models,  the  NSRL  and  NSCRL, 
which  consider  the  cases  of  a  subversive  element  and  a  political  campaign,  respectively. 
We  explore  the  behavior  under  a  wide  range  of  choices  for  parameters  and  initial 
conditions. 

In  Chapter  IV  we  study  a  stochastic  model  related  to  the  NSRL  model.  We  cite 
an  important  theoretical  result  which  shows  that  the  distribution  of  the  number  of 
individuals  in  each  class  of  the  NSRL  model  (and  all  models  similarly  posed)  when  the 
propagation  of  the  ideology  ceases  can  usefully  be  approximated  by  a  multivariate 
Normal  distribution  with  a  mean  coincident  with  that  determined  by  the  deterministic 
NSRL  model.  This  distributional  quality  is  conditional  on  fairly  weak  conditions  on  the 
initial  states.  This  makes  our  approach  a  compliment,  and  in  some  cases  an  alternative,  to 
agent-based  simulation  of  this  social  phenomena. 

The  most  important  result  from  our  work  is  not  the  NSRL  and  NSCRL  models 
themselves,  but  the  methodology  that  was  applied  in  creating  them.  The  framework 
developed  here  may  be  used  to  study  a  competitive  ideological  process  with  a  large 
variety  of  dynamics.  What  we  present  here  is  a  small  piece  of  what  is  possible. 

As  shown  in  Chapter  III,  Figure  10,  a  large  increase  in  initial  supporters  is 
required  to  overcome  a  relatively  modest  increase  in  contrarians.  This  shows  that  given 
the  rules  in  the  NSRL  model,  the  supporter  side  is  operating  at  a  significant  disadvantage. 
This  is  an  important  insight  into  Sea  Shaping  and  other  operations  that  the  US  is  currently 
involved  in. 

Secondly,  we  show  that  for  large  populations,  for  a  fixed  initial  number  of 

supporters,  parity  is  reached  if  less  than  10%  of  the  population  is  initially  counter.  In 

Chapter  III,  we  define  parity  as  the  ‘break  even’  point,  where  the  number  of  pro- 
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individuals  is  equal  to  the  number  of  counter-  individuals  at  the  end  of  propagation. 
Operationally,  the  counter-  side  will  always  need  at  least  10%  of  the  population  to 
prevent  spread  of  the  ideology  for  any  number  of  supporters. 


B,  FUTURE  RESEARCH 

In  our  opinion,  this  Thesis  contains  many  possibilities  for  both  applied 
mathematicians  and  probability  theorists.  Many  mathematical  subtleties  have  been  left 
unexplored.  Additionally,  in  introducing  a  new  class  of  models  we  have  painted  with  a 
small  brush  and  considered  only  a  few  special  cases.  Many  extensions  are  possible.  We 
feel  that  natural  continuations  of  this  work  would  include  the  following; 

1.  Application  of  Barbour’s  method  for  computing  the  variance  of  the  absorbing 
states  distribution. 

2.  We  have  only  considered  the  asymptotic  behaviors  of  the  models  proposed. 
We  have  not  developed  any  Measures  of  Effectiveness  (MOE’s)  which  take  the  time- 
dependant  behavior  of  the  system  into  consideration.  These  new  MOE’s  could  be  studied 
using  the  deterministic  models  proposed  in  Chapter  III. 

3.  Throughout,  we  have  assumed  constant  population  sizes  with  homogeneous 
mixing.  There  is  no  reason  that  the  theory  developed  here  could  not  be  extended  to  cases 
where  the  mixing  parameters  {f3-)  are  continuous  functions  in  time,  or  to  cases  with 

changing  population  size  due  to  individuals  joining  the  population,  leaving  the 
population,  or  both. 

4.  We  have  not  fully  considered  the  impact  of  mass  media.  We  propose  that  a 
mechanism  which  would  deal  with  the  effect  of  mass  media  would  be  to  use  an  artificial 
redistribution  of  the  population  at  some  time  in  the  model.  This  would  be  accomplished 
the  same  way  that  Eanchester’s  equations  with  reinforcements  are  modeled. 

5.  We  were  not  able  to  find  an  adequate  data  set  on  which  to  fit  parameters.  If 
one  could  be  found,  a  fitting  to  the  data  and  hypothesis  testing  would  be  an  important 
project. 
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APPENDIX  1:  IMPLEMENTATION  OF  ALGORITHIM  I 

(MATLAB) 


function  [ ] =Vectors (n, lo, Co,  plot) 

%%  Code  to  compute  the  joint  probability  distribution  of  the  absorbing 
%%  states  in  the  NSRL  model.  Developed  by  Roberto  Szechtman  with 
Harrison 
%%  Schramm. 

%%  Uses  equations  (10,  11)  from  Chapter  IV. 

%  Input  variables: 

%  n:  population  size  (integer) 

%  lo:  Initial  pro-spreaders  (integer) 

%  Co:  Inital  counter-spreaders  (integer) 

%  plot:  Logical  (real,  treated  as  boolean) .  A  value  of  zero  suppresses 
%  plots,  a  value  greater  than  zero  produces  plots. 

%  Output  Variables 

%  left  void  in  on  this  page,  but  any  of  the  following  could  be  placed 
in 

%  the  brackets  immediately  after  'function' 

%mean  s,  mean  r,  mean  c:  The  expected  value  of  these  at  the  end  state, 
%computed  directly  from  the  definition  of  expected  value 

%Var  S,  Var  R,  Var  C:  The  variance  of  these  at  the  end  state,  computed 
%directly  from  the  definition  of  varaince 

%  Joint:  The  joint  distribution  of  classes  R  and  C  at  the  absorbing 
state . 


s=n-Io-Co;  i=Io;  c=Co; 
pr=zeros (n+1 , n+1 , n+1 ) ; 
pr (s, i, c) =1;  %  Working  probabilities 

f in=zeros (n+1 , n+1 ) ;  %%  end  state  of  (S  >0  and  C>0)  [s,c] 

f in2=zeros (n+1 , n+1 ) ;  %%  (Temporary,  S  =  0)  (not  used  as  of  now) 

while ( sum ( sum ( sum (pr) )) >0 )  %%  while  any  probability  of  forward 
transition  exists. 

tmp=zeros (n+1 , n+1 , n+1 ) ;  %  last  state; 

[s, i, c] =ind2sub (size (pr) , find (pr>0) ) ;  %  provides  col  vectors  of 
S, i, c  prob  >0; 
r=n+l-s-i-c; 

total=(s>0) .*(i>0) .*s.*i+(i>l)*.5.*i.*(i- 
l)+(r>0) .*(i>0) .*i.*r+ (c>0) .*(i>0) .*i.*c; 
total= (total==0) +total; 

pl= (s>0) .* (i>0) . *s . *i . /total;  %vectors  that  contain  the 
probabilities  moving  forward  for  each  element. 
p2=(i>l)*.5.*i.*(i-l) ./total; 
p3=(r>0) .*(i>0) .*i.*r. /total; 
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p4=(c>0) .*(i>0) /total; 


sl=s(s>l);  il=i(s>l);  cl=c(s>l);  pi l=pl ( s>l ) ; 
if (size (si) >0) 

tmp (sub2ind (size (tmp) ,sl-l,il+l,cl) ) =tmp (sub2ind (size (tmp) , si- 
1, il+1, cl)  )+pll.*pr (sub2ind (size(pr)  ,sl,il,cl)  )  ; 
end 

sl=s(s==l);  il=i(s==l);  cl=c(s==l);  pi l=pl ( s==l ) ; 
if (size (si) >0) 

fin2 (sub2ind (size (fin2) , il+1, cl) ) =fin2 (sub2ind (size(fin2) , il+1, cl) )+pll 
. *pr (sub2ind (size(pr) ,sl,il,cl) ) ; 
end 

s2=s(i>2);  i2=i(i>2);  c2=c(i>2);  p22=p2 ( i>2 ) ; 
if (size (s2) >0) 

tmp (sub2ind (size (tmp) ,s2,i2-2,c2) ) =tmp (sub2ind (size (tmp) , s2, i2- 
2, c2) ) +p22 . *pr (sub2ind (size(pr)  ,s2,i2,c2)  )  ; 
end 

s2=s(i==2);  i2=i(i==2);  c2=c(i==2);  p22=p2 ( i==2 ) ; 
if (size (s2) >0) 

fin (sub2ind (size (fin) , s2, c2) ) =fin (sub2ind (size(fin) ,s2,c2) ) +p22 . *pr (sub 
2ind(size(pr) ,s2,i2,c2) ) ; 
end 

s3=s(i>l);  i3=i(i>l);  c3=c(i>l);  p33=p3 ( i>l ) ; 
if (size (s3) >0) 

tmp (sub2ind (size (tmp) ,s3,i3-l,c3) ) =tmp (sub2ind (size (tmp) , s3, i3- 
l,c3) )+p33.*pr (sub2ind (size(pr) ,s3,i3,c3) )  ; 
end 

s3=s(i==l);  i3=i(i==l);  c3=c(i==l);  p33=p3 ( i==l ) ; 
if (size (s3)  >0) 

fin (sub2ind (size(fin) ,s3,c3) )=fin (sub2ind (size(fin)  ,s3,c3)  )+p33.*pr(sub 
2ind(size(pr) ,s3,i3,c3)  )  ; 
end 

s4=s(i>l);  i4=i(i>l);  c4=c (i>l) ;p44=p4 (i>l) ; 
if (size (s4) >0) 

tmp (sub2ind (size (tmp) , s4, 14- 
1,  c4  +  l) ) =tmp (sub2ind (size (tmp) , s4, 14- 
l,c4+l) )+p44.*pr (sub2ind (size(pr) ,s4,i4,c4) ) ; 
end 


s4=s(i==l);  i4=i(i==l);  c4=c(i==l);  p44=p4 (i==l) ; 
if (size (s4) >0) 

fin (sub2ind (size (fin) , s4, c4+l) ) =fin (sub2ind (size(fin) ,s4,c4+l) )+p44. *pr 
(sub2ind (size(pr) ,s4,i4,c4) ) ; 
end 

pr=tmp; 

end 
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%now  take  of  the  mass  associated  with  s=0  and  i>0 
clear  pr; 

fin3=zeros (n+1, 1)  ; 

pr=fin2;  %%  fin2  holds  the  case  where  s  =  0,  i  >  0,  OO . 
while ( sum ( sum (pr ) ) >0 ) 

tmp=zeros (n+1 , n+1 ) ; 

[i, c] =ind2sub (size (pr) , find (pr>0)  )  ; 
r=n+l-i-c; 

total= (i>l)*.5.*i.*(i-l)+(r>0) .*(i>0) .*i.*r+(c>0) .*(i>0) .*i.*c; 
total= (total==0) +total; 

p2=(i>l)*.5.*i.*(i-l) ./total; 
p3=(r>0) .*(i>0) .*i.*r. /total; 
p4=(c>0) .*(i>0) .*i.*c. /total; 


i2=i(i>2);  c2=c(i>2);  p22=p2(i>2); 
if (size (12) >0) 

tmp (sub2ind (size (tmp) , 12-2, c2) ) =tmp (sub2ind (size (tmp) , 12- 
2, c2) ) +p22 . *pr (sub2ind (size(pr)  ,i2,c2)  )  ; 
end 

i2=i(i==2);  c2=c(i==2);  p22=p2 ( i==2 ) ; 
if (size (12) >0) 

fin3 (sub2ind (size(fin3) ,c2) )=fin3 (sub2ind (size(fin3)  ,c2)  )  +p22  .  *pr (sub2i 
nd(size(pr) ,i2,c2) )  ; 
end 

i3=i(i>l);  c3=c(i>l);  p33=p3(i>l); 
if (size (i3) >0) 

tmp (sub2ind (size (tmp) , 13-1, c3) ) =tmp (sub2ind (size (tmp) , i3- 
1, c3) ) +p33 . *pr (sub2ind (size(pr)  ,i3,c3)  )  ; 
end 

i3=i(i==l);  c3=c(i==l);  p33=p3 ( i==l )  ; 
if (size (i3) >0) 

fin3 (sub2ind (size(fin3) ,c3) )=fin3 (sub2ind (size(fin3) ,c3) )+p33. *pr (sub2i 
nd(size(pr) ,i3,c3) ) ; 
end 

i4=i(i>l);  c4=c(i>l) ;p44=p4 (i>l) ; 
if (size (14) >0) 

tmp (sub2ind (size (tmp) ,i4-l,c4+l) ) =tmp (sub2ind (size (tmp) , i4- 
l,c4  +  l) )+p44.*pr (sub2ind (size(pr)  ,i4,c4)  )  ; 
end 

i4=i(i==l);  c4=c(i==l);  p44=p4 (i==l)  ; 
if (size (i4) >0) 

fin3 (sub2ind (size(fin3) ,c4+l) )=fin3 (sub2ind (size (fin3) , c4+l) ) +p44 . *pr (s 
ub2ind (size(pr) ,i4,c4) ) ; 
end 

pr=tmp; 

end 

%%  fin3  holds  i  =  0; 
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fin;  fin3; 

mean  s=sum ( ( 1 : n+1 ) *f in) ; 

Var  S  =  sum ( ( 1 : n+1 ) . ^2  *  fin)  -  mean  s^2; 
mean  c=sum (fin*(l:n+l) ')+(!: n+l)*fin3; 

Var  C  =  sum (fin*((l:n+l).^2)')  +(l:n+l).^2*fin3  -  mean  c^2 
mean  r=  n  -  mean  s  -  mean  c; 
sum ( sum (fin) ) ; ,  sum ( sum ( f in3 ) ) ; 

%  fin  holds  s , c 

%%  fin3  holds  c  distribution. 

joint  =  zeros (n+1 , n+1 ) ;  %  r  and  c  matrix 

for  c=l : n 

for  r=l : n 

if (n+l-r-c>0) 

joint (r, c) =fin (n+l-r-c, c) ; %s>0,  r>0,  c>0 

end 

end 

end 

for  c=l : n 

j ointO (c) =f in (n+l-c, c) ;  %s>0,  r=0,  c>0 
joint (n+l-c,  c) =joint (n+l-c, c) +fin3 (c) ; %s=0,  r>0,  c>0 

end 

j oint00=f in3 (n+1 ) ; %s=0 ,  r=0,  c=n+l 

sum ( sum (joint) ) +sum (joint 0) +joint00; 

j  =  sum ( j  oint ' ) ; 

N2  =  ( (1 :n+l)  . ^2)  '  ; 

Var  R  =  (j*N2)  -  mean  r^2; 

if (plot  >0 ) 

handlel=  figure; 
mesh (joint) ; 
xlim ( [ 0  n+1 ] ) ; 
ylim ( [ 0  n+1 ] ) ; 

xlabel (' Number  in  class  C'); 
ylabel (' Number  in  class  R'); 
zlabel ( ' Probability '  )  ; 
end 
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APPENDIX  II:  MATLAB  CODE  TO  PERFORM  BINARY  SEARCH 
FOR  CRITICAL  PROPORTION  OF  CONTRARIANS 


function [num  rout]  =  estimate (rO, tol,  minN,  maxN,  Isize,  step) 

%  July  2006  Harrison  Schramm 

%  works  with  SIRC2.m  written  by  Harrison  Schramm  and  Carlos  Borges 

%  Determines  the  number  of  initial  contrarians  (ro)  to  achieve  parity 
%  for  a  wide  range  of  populations  and  a  fixed  lo 

%  Input  Parameters 

%  rO:  Initial  guess  for  the  proportion 
%  tol:  Numerical  tolerance  for  binary  search 
%  minN,  maxN :  Bounds  on  the  size  of  the  population 
%  Isize:  lo  (fixed) 

%  step:  Number  of  samples  to  compute  on  the  interval. 

%  Returns  rout,  which  is  a  vector  of  the  values  of  rc  found  for  the 
%  different  population  sizes  and 
%  num  which  is  the  vector  of  population  sizes 

num  =  linspace (minN,  maxN,  step) ; 

rout  =  ones ( 1 , length (num) ) ; 

maxruns  =  100;  %  prevents  infinite  loops 

for  i  =  1 : length (num) 

rlow  =  max ( 0 , rO- . 05 )  ; 
rhi  =  min(l,r0  +  .05); 
rout ( i )  =  rO ; 
ratio  =  2; 
goagain  =  true; 
count  =  0; 


while (goagain) 

rout(i)  =  (rhi  +  rlow)/2; 
z  =  num ( i ) ; 
y  =  rout  (i) ; 
m  =  SIRC2(y,z,  Isize); 

1  =  length (m) ; 
count  =  count  +1; 
ratio  =  m ( 1 , 4 ) /m ( 1 , 3 ) ; 
if (count  >maxruns) 
goagain  =  false 


end 

if  (abs (ratio  -  1)  <  tol) 
goagain  =  false; 


end 

if (ratio  <1 ) 

rlow  =  rout ( i ) ; 

end 

if (ratio  >  1 ) 
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rhi  =  rout  ( i )  ; 

end 


end 

end 


n  =  num; 

n; 

rout; 

handlel  =  figure () 
plot (n, rout) 

title (' Critical  proportion  of  contrarians  vs.  population  size') 

xlabel (' population  size') 

ylabel('ro,  proportion  of  contrarians') 


function  [y]  =  SIRC2 (r, n, io) 

O, 

o 

%2006  Harrison  Schramm  and  Carlos  Borges 

%  yd)  =  S 
%  y(2)  =  I 
%  y(3)  =  R 
%  y(4)  =  c 

%Initial  conditions 

IC  =  [((l-r)*n  -io) ; io; 0;  (r*n) ]  ; 

[t  y]  =  ode45 (@f,  [0; 100]  ,  IC)  ; 


q, 

o 

function  dydt  =  f(t,y) 

dydt  =  [-y (1) *y (2) 

y(2)*(y(l)-y(2)-y(3)-y(4)) 
y  (2)  *  (y  (2)  +y  (3)  ) 
y(2)  *y(4)  ]  ; 
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